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DETONATION  PROPERTIES  OF  PENTOLITE 

ABSTRACT 

The  hydrodynamic  theory  of  detonation  has  been  applied  to  determine 
the  detonation  velocity,  pressure,  and  product  composition  and  the 
isentropic  expansion  of  the  explosion  products  of  pentolite.  Hydrodynamic 
principles  were  applied  to  determine  the  initial  pressure  of  the  shock 
formed  in  air  after  completion  of  the  detonation  and  to  calculate  the 
pressure  and  velocity  distribution  behind  the  detonation  wave. 
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SYMBOLS 


dP 

n  =  sound  speed  = 

s 

D  =  detonation  velocity 

E  =  internal  energy  of  explosion  products  at  (T,p) 

E  =  internal  energy  of  explosion  products  at  300°  K  and  1  atm 

o 

E-E  =  internal  energy  of  explosion  products 

o 

p  =  pressure 


V* 

S 

T 

t 

u 

V 

7 

X 

t 

P 

a 


=  chemical  energy 

=  distance  from  point  of  initiation 
=  entropy 
=  temperature 
=  time 

=  mass  velocity 

=  specific  volume  (=  l/p) 

/(3£n  n\ 

b 

-  >^  /+ 

-  j./  z 


=  u/x 


=  density 

=  ^ 

I  pc 

^1 

=  u/c 


detonation  front  value  of  p 
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INTRODUCTION 


According  to  the  hydrodynamic  theory  of  detonation  the  properties 
of  the  detonation  front  in  an  explosive  are  determined  by  the  conditions 
of: 

(1)  conservation  of  mass,  momentum  and  energy 

(2)  the  Chapman-Jouguet  hypothesis 

(5)  the  equation  of  state  for  the  explosion  products 
(4)  the  distribution  of  the  chemical  elements  in  the  explosive 
among  possible  molecular  species  of  the  product  composition. 

Because  of  inadequate  knowledge  of  the  behavior  of  gases  and  solids  at 
the  high  temperatures  and  pressures  encountered  in  detonation  processes  the 
choice  of  a  suitable  equation  of  state  presents  a  formidable  problem. 

Several  equations  of  state  have  been  proposed  by  other  investigators  (for 
a  brief  summary,  see  Reference  1).  The  equation  of  state  used  here  is 
the  equation  used  extensively  by  Brinkley  and  Wilson  (Ref.  2).  With  their 
equation  of  state,  Brinkley  and  Wilson  obtained  reasonable  agreement  with 
observed  detonation  velocities  for  a  variety  of  explosives. 

Pentolite  has  been  used  extensively  in  the  experimental  study  of  blast 
parameters  (see  Ref.  5  for  compilation  of  blast  data)  and,  therefore,  it  is 
of  interest  to  determine  the  detonation  characteristics  of  this  explosive. 

The  explosive  50-50  pentolite  is  a  mixture  of  equal  weights  of  the 
explosives  PETN  and  TNT.  If  a  "molecular  weight"  of  250  grams  for  pentolite 
is  assumed  its  chemical  formula  may  be  written  symbolically  as  Q25  03_ 

q1|'  Additionally,  it  is  assumed  that  the  detonation  process  gives 
rise  to  the  reaction 

Cq^H^N^O^  =  n^CO+n^CO^  +  n^H^+n^OH  +  n^HgO  +  n^NH^+n^O  +  nQN2+n^02  (l) 


where  for  pentolite 


q  =  5.825  = 

r  =  5-910  =  2n^+n^^+2n^+^n^+hn^Q 
t  =  3.250  =  ng+n^+2nQ 
s  =  8.040  =  nj^+  2n2+nj^+n^+n^+2n^ 

The  distribution  of  the  atoms  among  these  molecular  species  (products) 
is  determined  by  assuming  several  chemical  reactions  and  the  conditions  of 
thermodynamic  equilibrium.  With  these  assumptions  and  the  above  atom-balance 
equations  a  solution,  consistent  with  the  conditions  imposed  by  the  equations 
describing  the  detonation  process^is  obtained. 

The  values  of  the  parameters  at  the  detonation  front  are  used  as  initial 
conditions  for  the  determination  of  the  iaentropic  expansion  of  the  explosion 
products.  The  data,  thus  obtained,  are  sufficient  forthe  determination  of 
the  initial  pressure  of  the  shock  wave  formed  in  the  neighborhood  of  the 
detonating  explosive.  The  results  are  sufficient  to  determine  a  progressive 
wave  solution  of  the  distribution  of  pressure  and  velocity  behind  the 
detonation  front. 

CONDITIONS  AT  THE  DETONATION  FRONT 

Across  the  detonation  front  the  usual  conditions  of  the  conservation 
of  mass,  momentum  and  energy  hold.  These  conditions  for  u^  =  0,  are:  (see 
Ref.  4) 


l2) 


(a)  p^D  =  Pj_(D-Uj^) 

(b)  p^+p^D^=Pj^+  p^(D-u^)^ 

(c)  E^-E^  .  A«.l/2(^  -  |^)(Pi+P„) 
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where  ^  represents  the  difference  between  the  heats  of  formation  of  the 
product  species  and  the  explosive,  i.e.,  ^=^n.H. 

f,  explosive 

The  are  given  in  Table  1  and  explosive  "taken  to  be  57 . 688kcal/250 
g  pentolite.  The  explosive  pentolite  in  its  condensed  state,  is  assumed 

•z 

to  be  at  rest,  of  density  I.65  and  in  air  at  one  atmosphere  pressure 

and  500°K.  In  these  Laboratories  the  normal  loading  density  of  pentolite 
is  1.65  g/cra^. 

Equations  (j)  are  to  be  supplemented  by  the  Chapman- Jouguet  hypothesis 
which  states  that  the  detonation  velocity,  D,  is  the  minimum  wave  velocity 
compatible  with  the  conditions  imposed  by  the  conservation  equations  and 
the  given  initial  state  of  the  explosive.  Equations  (3),  along  with  the 
equation  of  state,  detemine  a  one-parameter  family  of  solutions  and  the 
solution  for  which 


(^)  =  ;  where  v  =  l/p  C^) 

Hugoniot  o 

determines  the  conditions  at  the  wave  front  (Reference  4).  Equation  (4) 
is  equivalent  to  the  Chapman -Jouguet  hypothesis.  Furthermore,  it  follows 
that  D  =  u+c  for  a  Chapman-Jouguet  detonation  (Reference  4). 

EQUATION  OP  STATE 


The  equation  of  state  used  here  is  the  Brinkley-Wilson  equation 
(Reference  2)  viz., 

n  RT(l+xe^^)  n  RT  P(x) 

P(v-.1V^)  -  -S-js -  =  -S-5, -  (5) 

5^n  k 

where  x  =  - ^  —  (6) 

M(v-t]V  )‘ir 
S 


&  =  1/4,  p 

the  weight 
the  solids 


=  3/10,  n  is  the  number  of  moles  of  gaseous  products,  q  is 

fraction  of  the  solids  products,  v  is  the  specific  voliune  of 

s 

and  the  k^  are  covalurae  constants  and  are  presented  in  Table  2. 


11 


Assuming  that  M  =  25O  grams  eind  that  the  only  solid  product  in  the 
explosion  products  is  graphite  of  density  2.25  g/cm"^  and  that  the  volume 
change  of  the  graphite  is  neglible,  the  form  of  the  equation  of  state 

4 

becomes 


p(V  -  16/5  n, ,)  =  n  RTFCx)  ;  V  =  Mv 


(7) 


11'  g 

From  thermodynamic  considerations  and  equation  (7)  the  equation 
for  the  change  in  internal  energy  of  M  grams  of  explosion  products  is 


±u  r—  -1 

I"  \  S  (T)  -  E.°(500°K) 

i=l  L. 

+  l/4n  RTxe°*^^ 


4*  n 


r 


T 


J  (C) 

500 


p  graphite 

(8) 


dT 


where  E=°(T)  -  E.°(30Q°K)  and  '  C  °  dT  are  given  in  Table  3.  The 


300 

entropy  of  M  grams  of  the  explosion  products  is 


S/l 


.  =  ^  - ^  n  in  ~  -  n 

1  1=1  g 

-  lA  ^  +  I  toT 


82.051+  n 


in 


S. 


r 

-r 

10 


(9) 


where  (t)  is  given  in  Table  4  and  K  =  z:  n.K. .  Prom  equations 

•  ^  — T  ^  ^ 

-L— J_ 

(8)  and  (9)  the  Helmholtz  function  A  =  E  -  TS  becomes 


=  ^  '  n.A.  (t)  +  RT  y  n.  in  n. 

— *  11''  -5 — i  1  1 

1  1=1 

+  n^RT  |in  X  -  (l-e®'^^)  -  in(K/RT^^^^ 


^10) 
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CONDITIONS  FOR  THERMODYNfAMIC  EQUILIBRIUM 
The  general  conditions  for  thermodynamic  equilibrium  are: 


and 


=  0 

1'  1 

1 


(constant  T,V) 


where 


\ 

T,V 


(11 


(12) 


(15) 


and  where  the  v. 's  are  the  stoichiometric  coefficients  of  the  reaction 
1 

imder  consideration.  Now,  frCM  (lO)  and  (13) 

0.  =  A.°(t)  +  RT  £n  n.  +  RT 

-  n  RT  k^/K  +  RT  [in  x  (l-3°*^^)  -  in  (k/ET^^^)J  (14) 


.  n  ETP(x)  i  ) 


and  at  equilibrium 

i:  -  0  = 


.'.(A  “(T)  +  BTl  +  RpI^v,  in  n.-  n 
i'  i  '  '  '  ^  1  1  i 


K 


+  jin  X  -  (l-e^'^"")  -  in  K/ETp^'^^^jZlv^  +  n^  F(x)  i  Hv. 

Now  n\..°(T)  +  RtI  =  -  RT  in  K  where  K  is  the  ideal,  equilibrium 

constant  of  the  jth  reaction  and  is  dependent  on  the  temperature.  Thus, 
(15)  can  be  rewritten  as 

^v^in  n^  =  in  -*1^1:?"  p/ngP(x)  +  10/3(e®*^^  -1^ 


(15) 


0.3x  Z:V-»  p 

+  n  xe  ''  +  4= 


K  ■  RT  ^i^s^  ""i 
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If  t  .  is  defined  by  the  expression 


-  r*,  p  ^ 

-^n  5  .=  I 

J  L  n  F(x)  3 

D 


10, 


+n  xe 
g 


^  7..  5"  V  .k. 

0.5x^  11 


K 


+  P  Jv.v  m 
RT-^  1  ^ 


then,  from  (l6)  and  (l7) 


TT  (^)  ^  =  K  K.  (18) 

1  pj  j  J 

i 

jth  reaction 


which  is  the  desired  condition  on  the  number  of  moles,  n^,  in  the  jth 
reaction.  Values  of  were  obtained  from  the  equation  (Ref.  5) 


-in  K 


PJ 

F°-E 


AE  °  n  fr 

RT  i^ll  i- 

0  '' 


RT  J,\ 


•^i/ 


(jth  reaction) 


■i'he  values  of 


RT 


and  AE  were  obtained  from  references  (lO  and  11) 

°o 

and  the  calculated  values  of  K  .  are  given  in  Table  5*  In  'tke  determination 

P  J 

of  the  chemical  composition  of  the  explosion  product  mixture  the  following 
reactions  are  assumed 


j=(l);  COg-f  Hg  =  CO  HgO 
.j=(2)j  2C0  =  C  +  COg 

j=(3);  l/2No  +  =  NO 

C.  (-  d. 

j=(4);  2H2O  =  PHg-f  Og  (19) 

J=(5);  5H2+  CO  =  CHj^+  HgO 
J=(6);  I/2N2  +  5/2H2  = 

j=(7);  HgO  =  I/2H2  +  OH 


where  the  numerical  values  of  j  denotes  the  jth  reaction.  Hence,  from 
equations  (18)  and  the  reactions  (l9)  there  follows 


(17) 
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=  n^n^/n^n^ 
2 

IC.  =  n^n-, 

^  ^  X 

K.  =  n^n^/n^iiQ 


1/2 


IT  IT  n 
1).  3  ^9'  ‘"5 


/  o/^  N 


^*^5  '^l0^5/’^3  “l 

/  ^/2  l/2 

Kg  =  ng/n^-"  ng  ' 

Kr  = 


and  from  (20;  and  the  atom-halarice  equations  (2)  follows  the  system  of 
equations  for  the  determination  of  the  n^'s: 


-1+  ^  l+il-K^G 

2, 


where  G=2s-r+2n ,+3n^+hnT --n,  -2n„-4n. 

3  o  iu  4  /  i 


^2  ' 

“3  =  V5/''2^  ^ 


=  r-s+n^+2n2+n^+2n^-2nj-3ng-4n^l^, 


n^  =  n,  l/n^  , 


“7  = 


-7  1/  *^cr^’z 

:?  •  u  D  o 

- 


n. 


K,  , 


ng  =  l/2(t=ng-n^)  , 

"9  =  Vn/ 

“10  " 

nil  ^  ‘l‘''l'“2"'l0 


15 


(21) 


.h  - — ' 


.^j 


where  n,  refers  to  the  niamber  of  moles  of  CO,n  the  rraiuber  of  moles  of  COo^ 

JL  ^  d. 

etc=  (See  equation  2). 


The  conditions  at  the  detonation  front  are  determined  by  solving 
equations  (3) > (4) , (6) , (7) , (8) j  and  (2l).  To  obtain  the  conditions  at 
the  wave  front  the  Hugoniot  curve  in  the  neighborhood  of  the  detonation 
point  must  be  determined.  The  followir^  iterative  procedure  was  used 
to  determine  the  Hugoniot  curve.  Starting  with  any  value  of  T>3000°K 
and  appropriate  values  of  p^°^,  and  n^*s,  consistent  with 

the  atom-balance  equations,  the  equation 

11 


(AQ)^°I  n. '°'e.(T) 
_  _ 


=  0.49680  x'  'e 


V^6/3  n^^ 


^■22) 


f  1^  5 

obtained  from  equations  (3),(7)j  and  (8), is  solved  for  V'  '^cm  /25O  g  pentolite. 


T 


The  symbol  e.(T)  denotes  E  °  (T)-E.°  (500°K)  for  i  =  1,2,..,, 10  and  C 

(1)  ^  500  ^ 

for  i  =  11.  p  is  then  obtained  from  the  equation  of  state 


dl 


82,054  ng^°^T 


16/3 


(7) 


With  these  values  of  p^^^,  and  T,  x^°^and  n^''^'’s  improved  values. 


^  1) 

n.''  of  the  n.*s  are  obtained  by  solving  equations  (l8)  and  (2l)  and 

(1) 


X '  '  is  obtained  from  the  relation 


.(1) 


K 


(1) 


-  K°  i 

"11  J 


(6) 


( 1 )  ( 1 ) 

where  K'  =  ^  n.  '^k..  By  increasing  superscripts  by  one  in  equation  (ll) 

1  ^  ( 2^ 

and  improved  value  V'  '  is  obtained  from  equation  (ll)  and  improved  values 
of  the  remaining  quantities  are  obtained  as  above.  The  process  is  continued 

until  the  convergence  criteria  |  ^  lO”*^  |  lo"^ 

and  c  10 "^x'^'  and  j  n^'"J''-n.^  are  satisfied.  For 

these  calculations  cn  and  p  are  chosen  such  that  6  —  cn  ^  P  9- 
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The  converged  values  of  p,  V,  x,  and  T  determine  a  "point"  on  the 

Hugoniot  Giirve.  By  assuming  different  values  of  T  additional  points  on 

the  Hugoniot  can  be  determined  by  the  above  process.  The  detonation 

conditions  are  detexTuined  by  locating  the  pointy  p^  V  and  T  on  the  Hugoniot 

/dp.  P-Pn 

GXirve  which  satisfies  the  condition  that  ^  where  p  =1  atm, 

dV  H  V-V  o  ' 

5  ^ 

=  151 -5  cm''^/250  g  pentolite  and  where  the  differentiation  refers  to  the 

Hugoniot  cxirve.  The  computed  detonation  values  are 

p^  =  231250  atm 

-  116.19  cm^/250  g  pentolite 
T^  =  3567. 7°K 
D  =  7607  m/sec 

Gehring  and  Dewey  (Ref.  7),  report  experimentally  determined  values  of 
p^=  257000  atm 

and  D  =  7^20  m/sec 

for  pentolite  of  density  I.655  g/cm"^. 

ISENTROPIC  EXPANSION  OF  THE  EXPLOSION  PRODUCTS 


In  genereil,  the  explosion  products  are  a  mixture  of  gases  and  solids 
and  when  such  a  mixture  expands  Isentropically  the  chemical  composition 
will  change  due  to  varying  equilibrium  conditions.  In  determining  the 
isentropic  expansion  it  is  assumed  that  thermodynamic  equilibrium  exists 
throioghout  the  initial  phase  of  the  expansion. 


If  the  expansion  is  isentropic  then  the  condition  that 

J-C  £- 


(23) 


holds  throughout  the  expansion.  The  subscripts  m,  and  nn  indicates  that  the 
chemical  composition  of  the  mixt'ore  at  T^^  may  be  different  than  that  at  Tg* 

The  procedure  followed  in  the  calculation  of  the  isentropic  is  to  let  T  be 
the  independent  variable  and  to  calcilLate  the  pressxire,  volume  and  equilibrium 
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composition  at  T.  The  starting  conditions  for  the  calculation  of  the 
isentrope  are  the  detonation  front  values  of  T,  p,  and  V.  The  detonation 
front  values  are  indexed  by  subscript  one,  i.e.,  p^^,  Y^,  T^^,  x^,  n^  ^ 
i  =  1,2,,,., 11.  The  entropy  (T^)  of  the  mixture  of  explosion  products 

is  calculated  from  equation  (9)  and  Table  6.  It  is  assumed  that  this 
value  of  the  entropy  remains  constant  throughout  the  expansion.  If  Tp 
is  such  that  T,  -  100  T„  ^  T'  the  values  of  p,  V  and  the  n.  's  at  T^ 

i  —  i:  1.  X  ct 

can  be  determined  as  follows . 


An  approximate  value  of  the  equation  of  state  variable,  x,  at  Tg. 
is  obtained  from  equation  (23)  by  initially  assuming  that  the  chemical 
composition  remains  invariant.  With  this  initial  assumption  is 

obtained  by  solving  the  equation 


in  Xp^°^  +  10/3(1-0. 075Xp'°')  ^  =  (^  -  ^ 


*+  ip 


11 

-  IT  n.  , 
“g,i 


3  ^  1^ 

S^°(T^)  -  S^°(T2) 


1  +  in  X 


R 


and  then  Vp'"'  is  obtained  from  the  definition  of  x,i.e.. 


1 
(24) 


I"  n.k. 
/  \  1  1 

V 

2  To; 

T 


lb 


(o)  „  1/4  3  ""ll^l 

V  '  *  J ' 


Substitution  of  Vp'°^  and  Tp  into  equation  (j)  gives  With 

these  approximate  values  of  x,  p  and  V  at  Tp  the  ’.'alues  of  are  obtained 
from  equations  (I8)  for  each  of  the  reactions  (19)  and  then  from  equations 
(21)  the  corresponding  n^*s  are  obtained  by  Iteration.  From  equation  (9) 
and  the  definition  of  x  (6)  equation  (I5)  is  obtained. 


X  =  -  4.01324268  +  42  G 


(25) 


where  G 


,  US  11 

“i  r  <’'2)  -  r  U  “1  "1  +  ^ 

g  1=1  g  1=1 


y-16  n 
5"  ^ 

T^ 


1  ^  I  ; 

-n¥-  -  -  1*5240442558  = 

"  g  2 


n 
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The  above  calculated  n^’s 


(o) 


,  Pg'  and  Vg 


(o) 

(1) 


are  used  to  obtain  a  second 
/  1  \  /  ^\ 

,  Pg''  and  n^''‘^''». 


approximation  to  from  (25)  and  the  Vg'”',  ^2'~'  n^'~''»B  are  obtained 

as  above.  This  process  continues  until  convergence  of  the  values  of  the 
parameters  is  obtained.  The  values  calculated  at  Tg  thus  represent  a  "point" 
on  the  isentrope.  The  calcixLation  then  proceeds  from  Tg  to  T^  by  the  same 
process  as  Indicated  above.  The  expansion  with  varying  chemical  composition 
is  continued  until  T  =  l840°K  where  it  is  arbitrarily  assumed  that  for 
temperatures  less  than  l840°K  the  composition  remains  frozen.  The  results 
of  the  computation  are  presented  in  Tables  6  and  6a  and  Figure  1.  Table  6 
Includes  calculated  values  of 


1/2 


-V  - 


;  W  =/  p  dV  ;  a  =  f 

^1  ^1 


pc 


and 


E 


Eq  (eq.  4) 


CALCULATION  OP  THE  INITIAL  SHOCK  ERESSUEE  IN  AIR 

When  6in  explosive  is  detonated  the  detonation  wave  travels  through  the 
explosive  until  it  reaches  the  surface  of  the  charge.  Produced  at  the  surface 
is  a  shock  wave  in  air  and  a  rairef action  wave  in  the  explosion  products. 

The  shock  wave  moves  outwaird  into  the  air  and  the  rarefaction  wave  travels 
backward  through  the  explosion  products.  The  explosion  product  -  air 
boundary  moves  outwaird  also  but  at  a  slower  rate  than  the  shock  wave.  When 
the  wave  is  normal  to  the  surface  dynamical  considerations  require  that 
pressiire  and  mass  velocity  be  equal  across  this  boundary.  The  mass  velocity 
In  the  explosion  products  at  the  bovindary  is  given  by  (Ref.  4). 

(26) 

where  the  subscript  (l)  refers  to  the  detonation  front.  In  air,  the  velocity 
and  pressinre  are  related  by  the  Hugoniot  conditions;  i.e., 

u  =  u(p)  (27) 
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The  integral  in  equation  (26)  is  tabulated  in  Table  6  and  the  condition 
(27)  is  tabulated  by  Shear  and  Day  (Ref.  7)-  Using  these  data,  equations 
(26)  and  (27)  were  solved,  graphically,  and  an  initial  shock  pressure  of 
670  atm  was  obtained.  This  value  depends  on  the  temperature  where  the 
composition  is  frozen  since  when  the  composition  was  frozen  at  2590°K  the 
value  of  p  was  665  atm.  Sultanoff  (Ref.  8)  reports  an  experimental  value 
of  818  atm.  The  accuracy  of  this  measurement  is  uncertain  since  the  method 
by  which  it  was  determined  is  essentially  an  extrapolation  procedure  and 
is  dependent  on  the  method  of  data  reduction. 

Holt  (Ref.  12)  has  calculated  the  initial  pressure  of  the  shock  wave 
formed  in  air  at  the  surface  of  FETN.  He  reports  a  pressure  value  of  675 
atm.  Jones  and  Miller  (Ref.  13)  give  the  adiabatic  expansion  of  the  explosion 
of  TNT  at  a  loading  density  of  1.5s/cm'^.  The  data  reported  by  Jones  and  Miller 
are  sufficient  to  determine  the  initial  pressure  of  the  shock  formed  in  air 
at  the  surface  of  TITT.  The  present  author  made  approximate  calculation  from 
this  data  and  obtained  a  value  of  515  atm. 

PRESSURE  AND  VELOCITY  DISTRIBUTION  BEHIND  PLANE  AND  SPHERICAL  DETONATION  FRONTS 

The  specification  of  the  conditions  at  the  detonation  front  and  of  the 
pressure-volume  relation  during  the  isentropic  expansion  of  the  explosion 
products  provides  the  data  necessary  for  the  determination  of  the  pressure 
and  velocity  distribution  behind  the  front.  The  general  problem  was  considered 
first  by  G.  I.  Taylor  (Ref.  9)  who  applied  the  results  of  his  analysis  to 
obtain  the  pressure  and  velocity  profile  behind  the  detonation  wave  in  TNT. 
Taylor  assumed  that  the  detonation  velocity,  D,  is  constant  and  that  the  time 
necessary  to  establish  the  detonation  is  neglible.  With  these  assumptions 
all  parameters  behind  the  front  are  functions  of  X  =  r/t  where  r  is  the 
radial  distance  from  the  point  of  initiation  and  t  is  the  time  given  by 
R/D  =  R/(Uj^+c^)  where  R  is  the  distance  of  the  detonation  front  from  the 
point  of  initiation.  The  equations  describing  the  flow  behind  a  spherical 
front  are 
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4-  1,  4-  ^  -  ^P’^ 

- ~ 


Su  ,  „  3u  G  ^P 

3t  5?  =  -  T  ^ 


and  the  condition  that  p,  u,  and  p  depend  on  X  =  r/t  is  expressed  by 

L(u;  p;  p)  =  0  where  ^  =  |:£-  +  X  (29) 

Following  Taylor  and  introducing  the  variables 

5  =  u/Xi  t  =  u/c  and  Z  =  inX  (50) 

equations  (27)  and  (28)  become 


(1-0 


f  2C-(1.C)^^  f} 


.3^^-  (1=0^^ 
Hz  "  — r'.'2T2‘~^ 

U-U  r  -0 


where 


p  dc 

T  “dp" 
c 


Since  R  =  D  t  and  X  =  r/t;  there  follows  the  relation 

X/D  =  r/R  (52) 

Noting  that  Z  appears  only  as  a  differential,  dZ,  in  (20)  its  value  may  be 
taken  as 

Z  =  f n  r/R  (33) 

From  (32)  and  (53)  and  the  definitions  of  ^  and  I'  there  follows  the  relations 


21 


u/B  =  ^  x/R 
0  =  c/d  =  4  r/R 


(54) 


Since  c  has  been  determined  (Table  6)  0  may  be  taken  as  the  independent 
variable  by  using  the  trsmsformation  defined  by 


SO  ~  ~  ^  ^  ^ 

Thus,  there  follows  from  (31) 


dt  ^  2^-(l-^)t^f 

^  0(l.S)Vf 


and  where  Z  can  be  determined  from  (53)  and  (54)  or  from 


/  -T/-  \ 


(56) 


^  0(i-C)C 


(57) 


The  calculation  of  the  above  equations  proceeds  inward  from  the  front 
where  the  initial  values  of  Z,  't'  and  0  are: 


_  =  O.P^^T^ 

Pi  — 


1  =  0,50405 


01  =  c^/D  =  pjp^  =  0.76685 


7.  = 


n 


(58) 
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The  solution  of  (36)  is  tabulated  in  Table  7  and  P  and  u/D  are  plotted 
versus  r/R  in  Figure  2  and  3-  The  particle  velocity,  u,  becomes  zero  at 
r/R  =  0.43  and  the  mass  of  the  explosion  products  contained  between 
r/R  =0.43  and  r/R  =  1  is  approximately  233  grams  which  is  about  93  ^  of 
the  total  mass  of  explosion  products. 

For  comparison,  the  particle  velocity  and  pressure  distribution  vs. 
r/R  behind  a  plane  detonation  wave  are  shown  in  Figure  2  and  3*  The 
results  shown  for  the  plane  detonation  wave  are  those  which  correspond  to 
the  physical  model  of  a  detonation  wave  traveling  in  a  rigid  walled  tube 
closed  at  the  end  where  initiation  takes  place.  If  the  detonation  wave 
is  planar  then  behind  the  front  constant  values  of  u,  p  and  c  travel  with 
the  speed  u  +  c  and  u  is  given  by 


where  u^^  and  p^^  are,  respectively,  the  detonation  front  values  of  particle 
velocity  and  pressure.  The  integral  has  been  tabulated  in  Table  6  and  thus 
u  can  be  determined  for  various  values  of  p  and  the  corresponding  values 
of  p  and  c  are  known  from  the  adiabatic  expansion  (Table  6) .  Thus,  pressure 
and  particle  velocity  are  determined  as  functions  of  r/R  =  (u+c)/d  where 
r  and  R  are  defined  as  before. 

For  the  spherical  wave  front'  calculations  the  chemical  energy  released 
when  u  =  0  is  l.l62  kcal/g  and  the  temperature  is,  approximately,  2390°K. 

The  detonation  front  value  of  AQ  is  1.109kcal/g  and  when  the  composition  of 
the  explosion  products  is  frozen  AQ  is  1.186  kcal/g.  Since  p  and  u  are 
known  (from  36)  the  corresponding  values  of  E  -  can  be  determined  and 
the  total  energy,  carried  by  the  wave  motion  can.be  calculated.  Now 

^total”  ^kinetlc^  ^internal  where 


E 


Int 


1 

(E-E^)P(^) 


2 


o 
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E,  .  =  JwR- 
km 


■r 


pu^(r/R)^ 


1 

p( r/R)^d(r/R) 


From  the  solution  of  equation  (36)  and  the  properties  along  the  adiabatic 
the  calculated  values  of  E^^^  and  Ej^^^  are  1.063  kcal/g  and  0.084  kcal/g 
and  =  1.132  kcal/g  and  the  distribution  of  energy  between  internal 

and  kinetic  is  approximately  95  ^  internal  and  7^  kinetic. 

CONCLUSION 


The  calculated  values  of  the  detonation  pressure  and  velocity  are 
231550  atm  and  7807  m/sec  respectively.  These  values  agree  to  within 
2-1/2^  with  the  corresponding  experimental  values  reported  by  Gehring 
and  Dewey  (Ref.  6).  The  calculated  initial  pressure  of  the  shock  wave 


■f* n  >0  o  T  v' 

axi  0.0 


1 4*  T  ft  ^  1  r>.T  Tf-\  tn  4-V>^v-. 

a,u*Ju.o  “ .  . 


oiiaix  ocxxxv^x  x  o  c;;wpcx  a.xixc:xx  l 


value  of  818  atm. 


The  chemical  energy  released  at  the  detonation  front  is  1.109  kcal/g 
and  when  the  composition  is  frozen  the  chemical  energy  released  is  I.I86 
kcal/g.  The  total  energy,  carried  by  the  wave  motion  is  1.152  kcal/g 

and  its  distribution  between  internal  and  kinetic  energy  is  approximately 
93^  internal  and  7^  kinetic. 
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TABLE  1 


Heats  of  Formation,  (cal/mole) 


Chemical  Specie 

Index , i 

H. 

1 

Reference 

CO 

1 

26413 

9 

GOg 

0 

r\l.  nrs 

9 

«2 

3 

0 

9 

Q 

4 

/~\  -7  i~0 

-^P?o.y 

10 

^20 

5 

57802 

9 

NH, 

P 

6 

15469 

10 

NO 

7 

1 

01 

— < — lOOw 

9 

•RT 

"2 

8 

0 

9 

^2 

9 

0 

9 

CH4 

10 

17904 

T  n 
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TABU5 


Covolume  Constants  (Ref.  2) 


Species 


Index,  1 


k.cm^  °W^/Mole 
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TABLE  5 


Ideal  Internal  Energies  Relative  to  ,jO<J°K  (cal/inole) 


e.(T) 


Species 

CO 

COg 

“2 

OH 

H^O 

NH, 

J 

NO 

«2 

O2 

CH^ 

C 

■K- 

T°K 

*^1 

"3 

*^4 

•"‘6 

‘'8 

*^9 

^10 

^11 

300 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

iiOO 

499-67 

741.72 

495.82 

511.32 

611.11 

691.37 

514.70 

498.54 

511.46 

709.25 

245.82 

500 

1007.1 

1571-5 

995.60 

1018.3 

1242.0 

1455.3 

840.01 

1002.4 

1043.9 

1548-7 

564.45 

6oo 

T  rnO  r\ 

,  U 

f' 

3432.4 

Ih-97  ♦  0 

1  rnir  .if 

0  »  U 

1  Or\>->  ’t 
(  .  ( 

2580.3 

2 1  1.  r\ 

^ 

t  •  y 

2134.1 

Xp-LtJ  .  f 

1  ^r\r^  ). 

XUUU .  *+ 

A 

1566.2 

700 

2065.5 

2000 . 4 

2031.6 

3217.8 

2045.2 

2179.6 

3648.2 

800 

2620.7 

4441.3 

2507.2 

2545.5 

3291.5 

4193.7 

2708.6 

2589.6 

2778-5 

4896.1 

1826.1 

900 

5192.8 

5491.2 

3019.6 

3064.9 

4032.3 

5239.2 

3300.4 

3150.1 

5395.8 

6257.5 

2511.8 

1000 

5780.0 

6575-2 

3539.0 

3596.7 

4803.4 

6551.9 

3907.3 

3725.7 

4022.7 

7727-1 

2818.2 

1100 

4380.3 

7687.8 

4o66 . 6 

4136.1 

5604.8 

7515.4 

4438.3 

4515.1 

4662.9 

9290.4 

5540.9 

1200 

4993.1 

^24.9 

4603-4 

4685.4 

6435-9 

8752.5 

5159.3 

4916.8 

5312.6 

10950 

3877.2 

1300 

5616.2 

9982.7 

5150.5 

5249.3 

7295.6 

10072 

5797-8 

5529.5 

5970.6 

12660 

4424.3 

liiOO 

6248.0 

11157 

5707.9 

5820.9 

8182.5 

11435 

6445-1 

6151.5 

6655.8 

i4199 

4980.5 

1500 

6887.4 

12347 

6276.2 

64o4.7 

9094.9 

12842 

7099-8 

6781.6 

7307.8 

16291 

5544.1 

l600 

7533.4 

13550 

6854.8 

7000 . 1 

10031 

14273 

7761.0 

7419.2 

7986.0 

18188 

6114.2 

1700 

8185.2 

14764 

7443.2 

7600.7 

10989 

15731 

8425.6 

8065.2 

8670.2 

20124 

6689.3 

1800 

8842.6 

15988 

8040.8 

8212 . 3 

11968 

17192 

9091.7 

8712.8 

9560.2 

22102 

7267.9 

1900 

9504.2 

17220 

8648.2 

8835.7 

12965 

18683 

9769.4 

9566.8 

10056 

24117 

7849.4 

2000 

10170 

18460 

9264.8 

9461.1 

13979 

20187 

10445 

10U26 

10757 

26159 

8453.5 

2100 

io84o 

19707 

9889.5 

10100 

15010 

21729 

10707 

10689 

11464 

28231 

9020.8 

2200 

11512 

20960 

10522 

10742 

16054 

23248 

11805 

11556 

12177 

50326 

9611.2 

2300 

12188 

22218 

11165 

11389 

17113 

24798 

12491 

12026 

12895 

52445 

10204 

2400 

12666 

23432 

lioiO 

12044 

18184 

26359 

13179 

12698 

15618 

34580 

10798 

2500 

13546 

24751 

12465 

12704 

19266 

27942 

13870 

15574 

14547 

56752 

11594 

2600 

14229 

26023 

13126 

13372 

20358 

29512 

14558 

i4o52 

15081 

56906 

U994 

2700 

14914 

27500 

13795 

14039 

21461 

31092 

15255 

14752 

15820 

4lo86 

12595 

2800 

15601 

28580 

14466 

14715 

22572 

32666 

15952 

15414 

16564 

45284 

15198 

2900 

16289 

29663 

15145 

15391 

23692 

54225 

16649 

16097 

17314 

45494 

15802 

3000 

16979 

31150 

15830 

16075 

24820 

35807 

17345 

16785 

18068 

47713 

14408 

3100 

17670 

32440 

16520 

16762 

25958 

37457 

l8o46 

17470 

18826 

49946 

15016 

3200 

18363 

33733 

17215 

17455 

27097 

39268 

18750 

18159 

19590 

52181 

15626 

3300 

19058 

35028 

17916 

18155 

28278 

40954 

19447 

18849 

20357 

54431 

16237 

31-00 

19754 

36327 

18622 

18850 

29400 

42622 

20155 

19?4i 

21128 

56685 

16850 

3500 

2o45o 

37627 

19333 

19556 

30563 

44314 

20857 

20234 

21904 

56944 

17464 

3600 

21148 

38931 

20048 

20260 

31726 

46l66 

21568 

20928 

22684 

61210 

18080 

3700 

21848 

40236 

20768 

20969 

32899 

47695 

22272 

21623 

25466 

63487 

18696 

3800 

22547 

41544 

21494 

21685 

34072 

49407 

22985 

22520 

24252 

65764 

19514 

3900 

23248 

42854 

22223 

22401 

35254 

51113 

23693 

25017 

25045 

68052 

19952 

4ooo 

23951 

44166 

22956 

23126 

36436 

52834 

2440H 

23716 

25836 

70340 

20552 

lioo 

24654 

45480 

25694 

23840 

37624 

54559 

25114 

24415 

26652 

72631 

21172 

1200 

25558 

46796 

24436 

24576 

38816 

56276 

25851 

25115 

27450 

74951 

21795 

1300 

26065 

48il4 

25182 

25508 

40012 

58010 

26543 

25817 

28252 

77256 

224i4 

4400 

26769 

49434 

25933 

26028 

41212 

59699 

27260 

26519 

29036 

79530 

23056 

4500 

27476 

50756 

26687 

26771 

42414 

6i4i9 

27975 

27219 

29842 

8i845 

23658 

4600 

28184 

52080 

27445 

27509 

43620 

63190 

28698 

27926 

50650 

84159 

24281 

4700 

28692 

53405 

28208 

28252 

44826 

648^4 

29413 

28650 

51462 

86477 

24905 

48oo 

29601 

54735 

28973 

28991 

46o4l 

66656 

30159 

29556 

32274 

88785 

25529 

4900 

30511 

56062 

29",' 46 

29744 

47255 

68355 

30849 

50042 

53088 

91110 

26154 

5000 

51022 

57393 

50520 

30491 

48473 

70095 

31572 

50749 

53905 

95434 

26779 

T 

®11  =  ) 
3oo“k 
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TABLE  5a 

Ideal  Internal  Energy,  E*^  -  E^  cal/mole.  Relative  to  0"'K 


O 

O 

Ref 

CO 

1489.2 

10 

pn 

1658  s  6 

10 

^2 

i44o.5 

10 

OH 

1523.8 

11 

HgO 

1785.4 

10 

NH, 

5 

1823.8 

11 

NO 

1610.9 

11 

Ng 

1488.9 

10 

Or 

■i)im  c. 

1  A 

CH. 

1815.4 

11 

29 


TABI£  4 


Entropy  (ideal)  S^'^(T)/H 


Species 

CO 

CO2 

«2 

OH 

HgO 

NO 

"2 

°2 

CH, 

4 

C 

_o 

T  K 

S^(T) 

"5'  ' 

Si^(T) 

s^Ct) 

s 

b'“  ' 

S.j{T) 

Sy(T) 

Sjd’) 

s  . 

■10'"' 

a  fTl 
"11'"' 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

500 

23.782 

25.727 

15.726 

22.106 

22.722 

25.201 

25.555 

23.050 

24.686 

22.418 

0.697 

400 

24.792 

27.085 

16.752 

25.134 

25.894 

24.484 

26.585 

24.060 

25.714 

23.728 

1.053 

500 

25-585 

28.259 

17 . 516 

25.927 

24.825 

25.361 

27 . 197 

24.849 

26.554 

24.891 

1.409 

6oo 

26.245 

29.247 

18.158 

24.575 

25.069 

26.527 

27.873 

25.502 

27.227 

25.971 

1.753 

700 

26.816 

30.145 

18.703 

25.122 

26.292 

27.574 

28.459 

26.066 

27.830 

26.992 

2.082 

800 

27.525 

50.956 

19.177 

25.601 

26.903 

28.159 

28.979 

26 . 566 

28.366 

27.962 

2.392 

900 

27.779 

31 . 696 

19.598 

26.027 

27.460 

28.897 

29.447 

27.016 

28.849 

28.887 

2.679 

1000 

28.196 

32.376 

19.978 

26.414 

27.974 

29.585 

29.875 

27.426 

29.287 

29.771 

2  =  948 

1100 

28.571 

33.005 

20.527 

26.768 

28.454 

30.236 

30.227 

27.804 

29 . 690 

30.616 

3.197 

1200 

28.934 

53.589 

20.649 

27.095 

28.905 

30.867 

50.632 

28.154 

30.061 

31.424 

3.433 

1500 

29.265 

34.136 

20.949 

27.403 

29.331 

31.475 

50.969 

28.481 

30.406 

32.198 

3.654 

1400 

29.575 

34.649 

21.232 

27.690 

29.736 

32.053 

51.284 

28.787 

30.728 

32.958 

3.861 

1500 

29.866 

35.131 

21.498 

27.962 

30.121 

52.605 

51.580 

29.075 

51.031 

53.647 

4.055 

i6co 

30.140 

35.585 

21.750 

28.220 

30.490 

35.130 

51.860 

29.347 

31.316 

34.527 

4.242 

1700 

30.400 

56.016 

21.990 

28.465 

30.843 

53.631 

52.124 

29.604 

51.585 

54.979 

4.418 

1800 

30.646 

36.425 

22.219 

28.696 

51.181 

34.103 

52.572 

29.848 

51.840 

35.605 

4.584 

1900 

30.880 

56.815 

22.458 

28.920 

31.507 

54.565 

52.611 

30.080 

52.084 

56.207 

4.746 

2000 

51.105 

37 . 186 

22.649 

29.135 

31.820 

34.999 

32.836 

50.501 

32.516 

56.785 

4.896 

2100 

31.316 

37.541 

22.851 

29.339 

32.122 

35.427 

33.052 

30.513 

32.538 

57.545 

5.037 

2200 

31.520 

37.881 

23.046 

29.535 

32.413 

35.830 

53.258 

30.715 

52.752 

37.880 

5.173 

2500 

31.716 

38.207 

23.233 

29.725 

52.694 

36.222 

33.456 

50.910 

52.957 

38.398 

5.304 

2400 

31.904 

38.520 

23.415 

29.907 

32.966 

36.603 

33-646 

51.096 

35.154 

38.898 

5.435 

2500 

32.084 

38.822 

25.590 

30.084 

53.229 

56.962 

55.828 

51.276 

53.545 

39.581 

5.560 

2600 

32.258 

39.112 

23.760 

30.254 

33.484 

37.502 

54.003 

31.449 

33.529 

39.849 

5 . 661 

2700 

32.426 

39.392 

23.924 

30.419 

33.751 

37.654 

34.174 

51.616 

53.707 

40.301 

5.777 

2800 

32.588 

59.663 

24.084 

30.579 

33.971 

57.954 

34.337 

51.777 

35.880 

40.740 

5.887 

2900 

32.744 

39.925 

24.259 

30.754 

54.204 

38.258 

34.496 

51.933 

54.047 

41.165 

5.994 

5000 

32.896 

40.178 

24.389 

50.884 

34.450 

58.561 

34.648 

^.084 

34.210 

41.577 

6.104 

3100 

33.043 

40.423 

24.536 

31.051 

34.650 

58.857 

34.797 

52.250 

34.568 

41.979 

6.200 

5200 

35.185 

40 . 662 

24.679 

31.175 

34.864 

39.184 

34.941 

32.572 

54.521 

42.367. 

6.295 

3500 

33.324 

4o . 894 

24.618 

51.312 

35.072 

39.477 

35.079 

32.509 

54.671 

42.747 

6.391 

3400 

33.458 

41.118 

24.954 

31.446 

55.276 

39.753 

55.215 

52.643 

54.817 

43.115 

6.476 

3500 

33.589 

41.337 

25.087 

31.578 

55.474 

40.050 

55.347 

32.773 

34.959 

45.473 

6.577 

3600 

33.716 

41.550 

25.216 

51.706 

55-668 

40.501 

35.476 

52.900 

55.097 

43.825 

6.657 

3700 

33.839 

41.757 

25.543 

31.851 

35.856 

40.554 

35.600 

55.023 

55.235 

44.164 

6.748 

5800 

33.960 

41.959 

r-  1.  ^r-1 

31 . 954 

36.041 

*+0.014 

55.722 

55.145 

35.365 

44.497 

6.833 

3900 

34.078 

42.156 

25.589 

52.074 

36.221 

4l.o62 

55.841 

35.260 

35.494 

44.821 

6.914 

4ooo 

34.192 

42 . 349 

25.707 

52.191 

36.397 

4i.304 

35.957 

53. 374 

35.620 

45.158 

6.995 

4100 

34.505 

42.537 

25.824 

52.305 

56 . 5696 

41.544 

36.070 

55.4861 

55.7441 

45.448 

7.065 

4200 

34.414 

42.721 

25.937 

32.418 

56.7582 

41.771 

56.181 

53.5951 

35.8650 

45.751 

7.146 

!• 

54.521 

42 . 900 

n/' 

32.528 

36.9054 

42.003 

56.280 

33.7017 

35.9835 

46.047 

7.226 

4400 

54.626 

43.076 

26.159 

32.635 

37.0651 

42.219 

36.395 

55.8059 

36.0995 

46.336 

7.302 

4500 

54.728 

43.248 

26.267 

32.741 

37.2256 

42.436 

56.498 

35.9079 

36.2132 

46.619 

7.377 

4600 

54.828 

43.417 

26.373 

52.844 

57.5789 

42.652 

36.600 

34.0077 

56.3246 

46.898 

7.447 

4700 

34.927 

43.582 

26.477 

32.947 

37.5313 

42.853 

36.699 

34.1055 

36.4338 

47.170 

7.523 

}i  Pnr* 
‘T'J'-'V-' 

35.023 

X  '71.  "Z 

nil.  cnrn 
C.U. 

53.046 

57.6807 

56.797 

54.2013 

36.  54 10 

47.435 

7-598 

4900 

35.117 

43.902 

26. 680 

33.145 

37.8274 

43.261 

36.891 

54.2952 

36 . 6461 

47.697 

7.668 

50Q0 

35.209 

44.057 

26.779 

33.241 

57.9713 

43.453 

36.984 

54.3873 

36.7495 

47.954 

7.739 
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TABLE  5 


•jt 

Equilibrium  Constants 


T°K 

K 

K 

K 

K 

K 

K 

K 

I'l 

P2 

P3 

P4 

P5 

P6 

^7 

1000 

0 . 6946 

5.8009  (-1) 

7.4287  (-15) 

7.5515  (-21) 

3.9207  (-2) 

5.9808  (-4) 

5.0979  (-12) 

1100 

1.0096 

9.0169  (-2) 

3.0198  (-15) 

1.7215  (-18) 

3.2934  (-3) 

3.2391  (-4) 

1.1702  (-10) 

1200 

1.3674 

1.9505  (-2) 

6.6515  (-12) 

1.6030  (-16) 

4.1665  (-4) 

1.9465  (-4) 

1-5962  (-9) 

1300 

1.7548 

5.2861  (-3) 

9.1242  (-11) 

7.4785  (-15) 

7.2288  (-5) 

1.2582  (-4) 

1.4622  (-8) 

l400 

2.1567 

1.7544  (-3) 

0 

0  /  ■»  T  S 

T  /*\»'\  C  /  C  \ 

CD  ^  1 1  i  1  \ 

U  .  UdHU 

9.7'8g8  (-8) 

1500 

2.5719 

6.7776  (-4) 

6.0796  (-9) 

3.5633  (-12) 

4.3570  (-6) 

6.2218  (-5) 

5. 0815  (-7) 

l6oo 

2.9875 

2.9714  (-4) 

3.3630  (-8) 

4.3858  (-11) 

1.3948  (-6) 

4.6859  (-5) 

2.1515  (-6) 

1700 

.  3.3943 

1.4406  (-4) 

1.5230  (-7) 

4.0267  (-10) 

5.1155  (-7) 

3.6479  (-5) 

7.6855  (-6) 

1800 

3-7928 

7.6216  (-5) 

5.8290  (-7) 

2.8990  (-9) 

2.0974  (-7) 

2.9155  (-5) 

2^5862  (-5) 

1900 

4 . 1740 

4.3311  (-5) 

1.9556  (-6) 

1  ^  r\t~t  1  t  Ci\ 

i-.<yy(y  (-0) 

r\  ).£■/'■?  t  Ci  \ 

\-oj 

2.4052  (-5) 

6.5811  (-5) 

2000 

4.5405 

2.6088  (-5) 

5.7279  (-6) 

8.5484  (-8) 

4.6283  (-8) 

2.0088  (-5) 

1.6410  (-4) 

2100 

4.8871 

1.6551  (-5) 

1.5280  (-5) 

5.5555  (-7) 

2.4306  (-8) 

1.7164  (-5) 

3.7533  (-4) 

2200 

5.2166 

1.0995  (-5) 

3.7319  (-5) 

1.3155  (-6) 

1.3516  (-8) 

1.4906  (-5) 

7-9481  (-4) 

2500 

5.5292 

7-5837  (-6) 

8.3923  (-5) 

4.5506  (-6) 

7.9536  (-9) 

1.3112  (-5) 

1.5804  (-5) 

2400 

1-  n-,  /" f' 

5 .  OXDD 

r*  1,  /  /'"  \ 

?.4:)i;'j  {,-0) 

1  r-rr-tQ/\  /  1.  \ 

X. ( (by  v-h; 

*1  /  C  ' 

} .  000  c  f  /~i  S 

H.uuu:;  k  y 

rs/'r^r^  /  -7  \ 

2500 

6.0956 

3.9927  (-6) 

3.5299  (-4) 

3.5965  (-5) 

5.138?  (-9) 

1.0428  (-5) 

5.2861  (-3) 

2600 

6.3514 

2.9778  (-6) 

6.6648  (-4) 

9.1665  (-5) 

2.0789  (-9) 

9.3506  (-6) 

9.0150  (-3) 

2700 

6.5951 

2.3015  (-6) 

1.1998  (-3) 

2.1790  (-4) 

1.4280  (-9) 

8.4846  (-6) 

1.4798  (-2) 

2800 

6.8240 

1.8080  (-6) 

2.0697  (-3) 

4.8719  (-4) 

1.0066  (-9) 

7.7591  (-6) 

2.3423  (-2) 

2900 

7.0198 

1  r-/'\  /  /'N 

^-0) 

3-4453  (-3) 

1.0314  (-3) 

7.2806  (-10) 

7.1192  (-6) 

3.5940  (-2) 

5000 

7.2219 

1.1962  (-6) 

5.5292  (-3) 

2.0740  (-5) 

5.3994  (-10) 

6.6153  (-6) 

5.3515  (-2) 

5100 

7-3965 

9.8419  (-7) 

8.6464  (-3) 

4.0079  (-3) 

4.0697  (-10) 

6.1219  (-6) 

7.7998  (-2) 

5200 

7-5732 

8.2396  (-7) 

1.5081  (-2) 

7.3848  (-5) 

3.1574  (-10) 

5.7536  (-6) 

0.11033 

3500 

7-7111 

7-0315  (-7) 

1.9549  (-2) 

1.3206  (-2) 

2.4540  (-10) 

5.4503  (-6) 

0.155^ 

3400 

7.8711 

6.0012  (-7) 

2.7955  (-2) 

2.2705  (-2) 

1.9466  (-10) 

5.1301  (-6) 

0.20885 

3500 

7.9921 

5.2752  (.7) 

5.9598  (-2) 

3. 8041  (-2) 

1.5674  (-10) 

4.8792  (-6) 

0.28018 

3600 

8.1299 

4.5888  (.7) 

5.4879  ^2) 

6.1485  (-2) 

1.2849  (-10) 

4.6553  (-6) 

0.56859 

3700 

8.2363 

4.0680  (-7) 

7.4852  (-2) 

9.7638  (-2) 

1.0574  (-10) 

4.4245  (-6) 

0.47891 

3800 

8.5456 

3.6275  (-7) 

0.10040 

0.15063 

8.8352  (-11) 

4.2487  1-6) 

0.61252 

5900 

8.4211 

3.2420  (-7) 

0.13506 

0.22866 

7.4121  (-11) 

4.0706  1-6) 

0.77578 

4000 

8.5259 

2.9258  (-7) 

0.17279 

0.55759 

6.5299  (-11) 

5.9101  (-6) 

0.96766 

4100 

8.6069 

2.6410  (-7) 

0.22271 

0.49054 

5.4264  (-11) 

3.7728  (-6) 

1.1955 

4200 

8.6781 

2.4195  (-7) 

0.28288 

0.69984 

4.7029  (-11) 

5.6288  (-6) 

1.4618 

4300 

8.7585 

2.2218  (-7) 

0.35520 

0.98189 

4.0984  (-11) 

5.5282  (-6) 

1.7707 

4400 

8.8059 

2.0608  (-7) 

0.44c;  32 

1.3595 

5.5940  (-11) 

3.4io6  (-6) 

2.1293 

4500 

8.8488 

1.9185  (-7) 

0.54616 

1.8585 

3.1580  (-11) 

3.3070  (-6) 

2.5577 

4600 

8.8966 

1.7866  l-i) 

0.66617 

2.5008 

2.8074  (-11) 

3.2100  (-6) 

2.9991 

4700 

8.9420 

1.6707  (-7) 

0.80702 

3.5508 

2.4994  (-U) 

3.1142  (-6) 

3.5264 

4800 

8.9864 

1.5714  (-7) 

0.96849 

4.3674 

2.2481  (-11) 

3.0576  (-6) 

4.1103 

4900 

9.0122 

1.4854  (-7) 

1.1565 

5.6667 

2.0220  (-11) 

2.y683  (-6) 

4.7677 

5000 

9.0597 

1.4079  (-7) 

1.5670 

7.5179 

1.8314  (-11) 

2.8842  (-6) 

5.4935 

*See  (19)  of  text. 
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TABLE  b 


Therniodynaiiiic  Properties  of  the  Explosion  Products  During  Isentropic  , 
Expansion  of  the  Wave  Front  in  Pentolite  at  a  Jjoading  Density  of  1.65  g/cm'^ 


T°K 

P 

atiQ 

X 

r 

5 

p,g/cm 

Cjiii/sec 

-a,ni/  sec 

W 

g 

E  -  E 

0  g 

aq!^ 

g 

3567.7 

231347 

3.3663 

3.29 

2.1516 

f>«37 

0 

0 

1.505 

1.109 

3300 

211480 

3.2950 

3.34 

2.0942 

5845 

160 

0.068 

1.448 

1.120 

3200 

185228 

3.1399 

3.36 

2.0129 

5593 

387 

0.161 

1.365 

1.131 

3100 

161740 

3.0834 

3.34 

1.9331 

532? 

607 

0.246 

1.287 

1.158 

3000 

140602 

2.9746 

5.31 

1.8535 

5044 

826 

0.327 

1.211 

1.144 

2900 

121752 

2.8643 

3.26 

1.7740 

4761 

1040 

0.404 

1.138 

1.147 

2800 

104783 

2.7513 

3.20 

1.6935 

4482 

1255 

0.478 

1.067 

1.150 

2700 

89400 

2.6336 

3.14 

1.6108 

4201 

1472 

0.549 

0.999 

1.153 

2600 

75673 

2.5125 

3.07 

1.5266 

3926 

1690 

0.617 

0.934 

1.156 

2500 

63369 

2.3865 

3.00 

1.4399 

3655 

1912 

0.683 

0.870 

1.159 

2400 

52510 

2.2562 

2.90 

1.3509 

3382 

2136 

0.747 

0.809 

1.162 

2300 

42987 

2.1214 

2.82 

1.2596 

3120 

2363 

0.809 

0.750 

1.164 

2200 

54663 

1.9811 

2.72 

1.1654 

2863 

2596 

0.869 

0.693 

1.168 

2100 

27503 

1.8363 

2.62 

1.0686 

2612 

2833 

0.927 

0.639 

1.172 

2000 

21397 

1.6862 

2.50 

0.9688 

2367 

3077 

0.984 

0.587 

1.176 

1900 

16272 

1.5312 

2.38 

0.8661 

2131 

3328 

1,038 

0.538 

1.182 

1800 

12164 

1 . 57'+2 

2.36 

0.7657 

1969 

3555 

1.090 

0.490 

1.186 

1700 

9024.6 

1.2202 

2.25 

0.6728 

1751 

3818 

1.136 

0.444 

1.186 

1600 

6530.2 

1.0655 

2.15 

0.5808 

1565 

4061 

1.180 

0.400 

1.186 

1500 

4600.6 

0.9125 

2.05 

0.4911 

1390 

4308 

1.222 

0.559 

1.186 

1400 

3151.1 

0.7642 

1.92 

0.4057 

1230 

4558 

1.261 

0.319 

1.186 

1300 

2090.3 

0.6229 

1.81 

0.3256 

1086 

4812 

1.299 

0.282 

1.186 

1200 

1341.0 

0.4924 

1.71 

0.2530 

957 

5070 

1.334 

0.246 

1.186 

1100 

828.4 

0.3751 

1.61 

0.1891 

845 

5551 

1.568 

0.212 

1.186 

1000 

495.3 

0.2760 

1.52 

0.1361 

750 

5595 

1.400 

0.180 

1.186 

900 

283.1 

0.1940 

1.45 

0.09557 

667 

5859 

1.450 

0.150 

1.186 

800 

154.3 

0.1301 

1.39 

0. 06088 

598 

6129 

1.459 

0.121 

1.186 

700 

79.35 

0.0827 

1.35 

0.03747 

559 

64o4 

1.486 

0.0942 

1.186 

600 

57.95  0.0496 

1.33 

0.02101 

486 

6686 

1.512 

0.0685 

1.186 

500 

16.46 

0.0276 

1.32 

0.01150 

458 

6978 

1.556 

0.0442 

1.186 

400 

6205  0.0139 

1.32 

0.00549 

389 

7285 

1.559 

0.0214 

1.186 

300 

1875  0.0061 

1.52 

0.00233 

359 

7611 

1.581 

0.0003 

1.186 

TABUS  6a 


Temx)eratui‘e ,  Volume  and  Number  of  Moles/250g  Pentolite  of  Each  Explosion  Product 
During  Isentropic  Expansion  of  the  Wave  Front  in  Pentolite  at  a  Loading  Density  of  1.65  g/cnr 


T°K 

„  cm 

^ 2 50gPentol Ite 

CO 

COg 

Hg 

OH 

HgO 

NH, 

3 

NO 

«2 

°2 

CH, 

4 

c 

5567.7 

116.19 

3.509 

0.843 

0.000 

0.009 

2.946 

0.005 

0.089 

1.569 

0.000 

0.000 

1.673 

5300 

119.38 

3.285 

0.872 

0.007 

2.946 

0.003 

0.058 

1.584 

- 

- 

1.668 

5200 

124.20 

3.255 

0.901 

- 

0.005 

2.946 

0.004 

0.051 

1.598 

- 

- 

1.669 

5100 

129.52 

3.218 

0.928 

- 

0.004 

2.946 

0.005 

0.016 

1.604 

- 

- 

1.678 

5000 

134.88 

3.180 

0.952 

0.002 

2.945 

0.006 

0.008 

1.608 

- 

• 

1.693 

2900 

140.92 

3.154 

0.969 

0.001 

0.002 

2.945 

0.007 

0.004 

1.610 

- 

- 

1.702 

2800 

147.62 

3.126 

0.985 

0.002 

0.001 

2.941 

0.008 

0.002 

1.610 

- 

- 

1.714 

2700 

155.20 

3.094 

1.003 

0.003 

0.001 

2.958 

0.010 

- 

1.610 

- 

1.727 

2600 

165.76 

3.065 

1.020 

0.004 

- 

2.955 

0.011 

1.609 

- 

1.739 

2500 

173.63 

3.024 

1.044 

0.007 

- 

2.927 

0.015 

- 

1.608 

- 

- 

1-756 

2400 

185.06 

2.989 

1.066 

0.012 

- 

2.918 

0.015 

- 

1.607 

- 

- 

1.769 

2500 

198.47 

2.951 

1.091 

0.019 

- 

2.906 

0.018 

- 

1.606 

0.001 

1.781 

2200 

214.51 

2.903 

1.124 

0.031 

- 

2.888 

0.020 

- 

1.605 

- 

0.005 

1.795 

2100 

233.95 

2.844 

1.167 

0.048 

2.862 

0.022 

1.604 

0.006 

1.808 

2000 

258.05 

2.771 

1.223 

0.074 

- 

2.822 

0.024 

1.603 

0.011 

1.819 

1900 

288.64 

2.681 

1.297 

0.109 

- 

2.765 

0.025 

- 

1.602 

* 

0.021 

1.826 

1800 

326.50 

2.618 

1.351 

0.136 

- 

2.720 

0.025 

- 

1.602 

- 

0.050 

1.825 
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TABLE  7 


Flow  Properties  Behind  a  Spherical  Detonation  Wave  in  Pentolite 


At  a 

Loading 

Density 

of  1.65 

A 

g/cm 

r/R 

\jf 

u/D 

c/D 

P  g/  era 

p  . 
apm 

1.00000 

.235 

.304 

.253 

.767 

2.152 

231350 

.99985 

.2280 

.2990 

.228 

.762 

2.156 

225970 

.2139 

.2852 

.213 

.748 

2.094 

2ii480 

.9955 

.2010 

.2724 

.200 

.733 

2.053 

197960 

.9865 

.1892 

.2604 

.187 

.716 

2.013 

185250 

.9767 

.1782 

.2488 

.174 

.700 

1.973 

175140 

■  .9644 

.1679 

.2374 

.162 

.682 

1.933 

161740 

.9492 

.1581 

.2259 

.150 

.  664 

1.893 

150880 

.9515 

.1488 

.2144 

.138 

.646 

1.855 

l4o6oo 

.9104 

.1597 

.2026 

.127 

.628 

1.814 

130940 

.8868 

.1510 

.1905 

.116 

.610 

1.774 

121750 

.8609 

1  ooc 

.1781 

.105 

.592 

1.734 

113030 

.8318 

.1159 

.  1650 

.095 

.574 

1.694 

104780 

.7992 

.1049 

.1507 

.084 

.556 

1.652 

96879 

.7659 

.0955 

.1356 

.073 

.558 

1.611 

89400 

.7255 

.0856 

.1193 

.062 

.520 

1.569 

82341 

.6841 

.0748 

.1017 

.051 

.505 

1.526 

75673 

.6592 

.0626 

.0824 

.o4o 

.486 

1.483 

69325 

.5888 

.0482 

.0606 

.028 

.468 

1.440 

65370 

.5297 

.0500 

.0555 

.016 

.450 

1.396 

57778 

.4550 

.0058 

.0061 

.003 

.455 

1.551 

52510 

.4500 

.0000 

.0000 

.000 

.450 

T  XliO 

u.  • 

51499 
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TABLE  8 


Flow  Properties  Behind  a  Plane  Detonation  Wave  in  Itentolite 
At  a  Loading  Density  of  I-65  g/cm'' 


D  =  7607  m/sec 


r/p. 

u/d 

c/d 

^atni 

T'K 

1.0000 

.23315 

.76685 

231350 

3367.7 

.96096 

.21260 

.74856 

211480 

3306 

.93071 

,19797 

.75274 

197960 

3250 

.90004 

.18560 

.71644 

185230 

5200 

.86890 

.16956 

.69954 

173140 

5150 

.83740 

.15551 

.6^09 

161740 

3100 

.80545 

.14123 

.66420 

150880 

3050 

.77349 

.12737 

.64612 

l4o6oo 

3000 

.74147 

.11362 

.^785 

150940 

2950 

.70965 

.09987 

.60978 

121750 

2900 

.67813 

.086li 

.59262 

115050 

2850 

.64908 

.07495 

.57415 

104780 

2800 

.61449 

.05849 

.55600 

96878 

2750 

. 58265 

.04456 

.53809 

89400 

2700 

.55092 

.03062 

.52050 

82541 

2650 

.51942 

.01661 

. 50281 

75675 

2600 

.48797 

.00242 

.48555 

69323 

2550 

.48262 

.00000 

.48262 

68280 

2541.5 
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FIGURE  2 
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